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In order to detect the gravitomagnetic clock effect by means of two counter-orbiting 
satellites placed on identical equatorial and circular orbits around the Earth with radius 
7000 km their radial and azimuthal positions must be known with an accuracy of 5r = 
10~^ mm and Srp = 10~^ mas per revolution. In this work we investigate if the radial 
and azimuthal perturbations induced by the dynamical and static parts of the Earth' s 
gravitational field meet these requirements. While the radial direction is affected only 
by harmonic perturbations with periods up to some tens of days, the azimuthal location 
is perturbed by a secular drift and very long period effects. It results that the present 
level of accuracy in the knowledge both of the Earth solid and ocean tides, and of the 
static part of the geopotential does not allow an easy detection of the gravitomagnetic 
clock effect at least by using short arcs only. 

1. Introduction 

Latest years have seen great efforts, both from a theoretical and an experimental 
point of view, devoted to the measurement of the general relativistic Lense-Thirring 
effect^^^ in the weak gravitational field of the Earth by means of artificial satellites. 

Among the satellite-based experiments recently proposed, one of the most in- 
teresting is devoted to the detection of the the gravitomagnetic clock effect'^^'*'^. 
It consists in the fact that two clocks moving along pro- and retrograde circular 
equatorial orbits, respectively, about the Earth exhibit a difference in their proper 
times which, if calculated after some fixed angular interval, say 2tt, amounts to: 

(t+ - r^h=2. ^ 4^ ^ 10-' s, (1) 

where and are the rotational angular momentum and the mass, respectively, 
of the Earth; c is the speed of light in vacuo. In^'^ it has been shown that for an 
orbit radius of 7000 km the radial and azimuthal locations of the satellites must be 
known at a level of accuracy of Sr < 10^^ mm and Scj) < 10~^ mas per revolution. 
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However, the studies conducted up to now on the feasibihty and the error budget 
of such an experiment are stiU prehminary"''®. 

In this paper we shaU investigate in a quantitative manner the systematic errors 
induced on the radial and azimuthal locations by the Earth solid and ocean tidal 
perturbations and by the static part of the geopotential. The paper is organized 
as follows. In Section 2 and 3 the radial and azimuthal perturbations, respectively, 
induced by the most relevant tidal constituents are investigated. In Section 4 the 
radial and azimuthal perturbations generated by the static part of the geopotential 
are worked out. Section 5 is devoted to the conclusions. 

2. The radial error induced by the Earth solid and ocean tides 

According to^'^ , the position perturbations in the radial direction can be expressed 
in general as: 



Ar = ^ Aa^ + i [ (eAa + aAe)^ + {aeAM)% (2) 

where a, e and A4 denote the satellite' s semimajor axis, eccentricity and mean 
anomaly. In eq.(^ the perturbation amplitudes are the rss values of the perturba- 
tions and small eccentricity approximations have been extensively applied. 

Since the difference in the proper orbital periods to be investigated is integrated 
over 2iT with respect to the azimuthal angle 0, as viewed by an inertial observer 
fixed with the distant quasars, we shall consider only the long period perturbations 
averaged over an orbital revolution. This is accomplished by assuming those values 
for the indices l,p, q^'^'^ which satisfy the relation I — 2p + q = 0. Since the tidal 
perturbations on the semimajor axis a are proportional just to I — 2p + q, all the 
terms in (Aa)^ and Aa of eq.(0) vanish and it reduces to: 



Ar = ^ [ (aAe)2 + (aeAM)^]. (3) 
For e = eq.(^) becomes: 

Articles = a Aetides- (4) 

For a constituent characterized by a given set of indices {I, m, p, q}, the first order 
tidal perturbation amplitude for the eccentricity turns out to bef 

Ae,^,,o.-^-^^^P^f^, (5) 

where Fimp and Gipq are the inclination and the eccentricity functions, respectively, 
as can be found in^^. 



*This holds for I — 2p + q = 0. Note that for the eccentricity there are no second order, indirect 
perturbations due to the oblateness of the Earth^'^"'^'^ , contrary to the node, the perigee and the 
mean anomaly. 
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Eq. (|^) allows to obtain a preliminary insight into those perturbations which, for 
a given set of indices I, m, p, q, vanish. Note that, for e = 0, eq.(||) could become 
singular. Concerning this problem it must be considered that, since Gipq oc g' ^ 
the behaviour of c>c e' "^'^^ is crucial. If \ q\ — 1 = k > 0, i. c. q > 1 or q < — 1, 



then for a circular orbit = and the perturbation vanishes. If | q| — 1 = 0, 
i.e. q = ±1, then — const. Problems may arise only if g = 0, but we shall see 
that, in general, in the cases in which q takes such value, I — 2p — also holds so 
that the perturbations identically vanish with no regards to the eccentricity or the 
inclination of the satellite. 

Let us start with the tidal perturbations of even degree. For I = 2n, p = 0, ...,l 
and I — 2p + q = 0,the allowable values for q satisfy the above stated conditions so 
that we can conclude that there are no radial tidal perturbations of even degree. 
Since for the solid Earth tides we consider only the / = 2 constituents, this result 
rules out their possible influence on the radial error budget in the gravitomagnetic 
clock experiment. 

Now we shall consider the odd degree case. For I — 3 there are no problems 
because q never vanishes. Moreover, for p = 1, q = —1 and p = 2, q = 1 — ^ = 1 
and l — 2p = ±1 while for the other sets of indices the perturbations vanish because 
— 0. When we consider the inclination functions corresponding to the indices 

G 

for which — ^ and I — 2p differ from zero, i. e. i^smi and i^3m2, and evaluate them 
for i = we find that only ^311(1 = 0) = — |. 

So we can conclude that the radial direction is perturbed only by the ^ = 3, m = 
1, p=l, 1 ocean tides. 

The full expression for the eccentricity perturbation amplitude due to ocean 
tides^'^'', in our case, is given byf 



AirGp^ f Re, 1+2 A + h ^ , - 2p)Fi„,pGipq 
^eimpqf = :+ { — ) (7^7—7) Gi„,f [ J, (bj 



where: 

• G — 6.67259 • 10^^ cnv" g^^ s^^ is the Newtonian gravitational constant"'^^ 

• n = is the satellite mean motion; GM® = 3.986 • lO^" cm^s'^ " 

• Pw — 1.025 g cm~^ is the water density^"* 

• i?0 = 6378 • 10^ cm is the Earth' s equatorial radius 

• k[ is the load Love number. 4 = -0.3075, fcg = -0.195015^1^ 

• C^^p cm is the ocean tidal height as can be found in EGM96 model^^ 
''yinpqf = {l-2p)uj+{l-2p + q)M+mtl+{j2-m) s + ja h + jip + j^ iV'+jg Ps in 
which uj and are the satellite' s perigee and node and the integers ji, i — 2, ..6 refer 
to the Doodson number^^'i^ by which each tidal constituent is classified. For the 
astronomical longitudes s,...ps see, e.g.^^. RecaU that j^^p^f has to be evaluated 
on the chosen reference orbit. 



tHere we shall consider only the prograde waves^'^. 
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For / = 3, TO = 1, p=l, q = -1, and by putting j^^p^f = eq.(^) becomes: 

Ar3ii„i/ = (8.80 • 10^5 x a'^^^ x Ppew x (7) 

Among the tesseral tides, the Ki (165.555) is by far the most important in perturb- 
ing the near Earth satelhtes' orbits^"'^^. So it seems reasonable to calculate eq.(^ 
for it in order to obtain an upper bound in the order of magnitude of the tidally 
induced perturbations on Ar. For such a constituent: 

C^^{Ki) = 0.95 cm, (8) 

and 73ii_i(-f^i) = ui + If we assume as reference orbit a secularly precessing 
ellipse "'^^j we obtain: 

Pp,,, = T-^= a'/' = (4-7639 • 10-25 cm-y^s) x a^^ (9) 

Eqs.(^-(^ tell us the important feature that for Ki, 1 — 3, m — 1, p ~ 1, q = —1 
the perturbation amplitude is independent from the satellite' s semimajor axis. For 
a = 7000 km and C2,o = —J2 = —0.00108261 we obtain Ppert = 50 days. By using 
eqs.(^)-(^) in eq.(|^) we obtain: 

Ar3ii_i(i^i) = 39.849 cm. (10) 

According to EGM96 model^^, the percentage error on C^(i^i) amounts to 5.2%; 
this yields 6r3ii-i{Ki) ~ 2.07 cm. Despite the amplitude of this long period mis- 
modeled perturbation is 2 orders of magnitude greater than the maximum allowable 
error Sr„iax = 1 • 10^^ mm, it must be noted that its period Ppert amounts to only 
50 days. This implies that if an observational time span Tabs which is an integer 
multiple of Ppert, i- e. some months, is adopted the tidal perturbative action of Ki 
can be averaged out. 

3. The azimuthal error induced by the Earth solid and ocean tides 

Concerning the angular variable which defines the position of the satellite on the 
orbit, for an equatorial, circular orbit it seems reasonable to adopt for its rate of 
change: 

(j) = uj + ncosi + M. (11) 

In it -I- f2 cos i represents an angular velocity around the direction of the orbital 
angular momentum^^; it is valid for any inclination angle i. In order to account 
for the fact that the orbit is circular we add to it Ai. See also^'^. About the 
perturbations on the latter Keplerian orbital element, it turns out that'^^ in AA4 
one has to consider also the indirect perturbations on the mean motion n due to the 
cross coupling with the semimajor axis a. Since they are proportional to I — 2p + q, 
they vanish when only long period perturbations are considered, as is the case here. 
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The perturbation amplitudes on the node, the perigee and the mean anomaly are 
proportional to! 

AQlmpqf OC , : Glpq , (12) 

Vl — e^smi ai 
Vl - dGipq cosi dFimp 

i^^lmpqf OC ; rirnp , (j/pq 7: , (ioj 

e de VI - sm i a« 

AMhnpqf CX ^^imp + 2(^ + 1) Flr„p Glpq. (14) 



By assuming 1 — = Vl — e-^ ~ 1 for e — > 0, with the aid of eqs.(|Tl|)-(|l4|) we 
obtain: 

A4>lmpqf OC 2{l + \)Fimp Glpq. (15) 

As already done in the previous section, cq.(p^ can be used in order to forecast 
which perturbations will vanish. 

For / — 2 only the combination Z = 2, p=l, g = yields a nonzero eccentricity 
function: G210 = (1 ~ e^)^'^/^ — 1. Among the corresponding inclination functions 
^2mi, for z = we have F201 = —1/2. The same conclusion holds also for / = 4 
with G42o(e = 0) = 1 and i^402(i = 0) = 3/8. So we can conclude that for I = 2,4 
only the zonal tides, both solid and ocean, cause nonvanishing perturbations on the 
satellite' s azimuthal variable. 

Concerning the odd degree perturbations, they all vanish since for ^ = 3, 5, g 
is always nonzero, so that, since Gipq oc e' 'I, for circular orbits all the eccentricity 
functions vanish. The conclusion is that the odd part of the ocean tidal spectrum 
does not induce systematic errors on the satellite's azimuthal variable. 

For a given constituent of degree order m and frequency / the full expressions 
for the solid Earth and ocean tidal perturbation amplitude (progressive waves only) 
are, respectively: 

^<l>tiipqf = -Y^ ( — Aim kl°2f Hr [2{l + 1) Fir^p Glpq], (16) 

^c^ZpTf - ( — (^) CLf + 1) Fi^p Glpq], (17) 

where: 

• g = 978.0327 cm is the acceleration of gravity at the surface of the Earth as 
if it was perfectly spherical^^ 

• ilmpqf ~ ilmpqf 

4 _ I2I + I (i~m)! 

• ^Im — Y 47r (l+m)\ 

• k[2^ is the Love number for the free space potential^"''^"' 

• are the Doodson coefficients with a different normalization^^^^^ 



■l-The second order, indirect perturbations will not be considered here since it can be demonstrated 
that they vanish in this case. 
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We shall start by considering the three most relevant I = 2 m = zonal tides: 

• 18.6 - year (055.565); Ppert = 6798.38 days; = 0.315; = 2.792 cm 

• 9.3 - year (055.575); Ppert = 3399.19 days; fc^o' = 0.313; = 2.72 • 10"^ cm 

• Sa (056.554); Ppert = 365.27 days; = 0.307; = -4.92 • IQ-^ cm; k'2 = 
-0.3075; C+^f = 2.54 cm 

For 1 = 2 eqs.(|6l)-(0) becomes: 

^^ocean ^ (_4.707 • lO" Cm^/^s^l) X fl-^/^ X Ppert X C,+ ^. (19) 

Note that, since for the I = 2 zonal tides Ppert does not depend on the satellite' 
s semimajor axis but only on the astronomical arguments, lS.(j)impqf depends on 
the orbit' s radius through a^^/^, contrary to A r{Ki), as shown in the previous 
Section. For a ~ 7000 km we have: 

• A0(18.6 - year) = -4.431 • 10^ mas 

• A(^(9.3 - year) = -214.4 mas 

• A(j){Sa) = 408 mas (solid); 857.6 mas (oceanic) 

The zonal tidal perturbations on the satellite' s azimuthal location are partic- 
ularly insidious not only because their nominal amplitudes are up to 6 orders of 
magnitude greater than the maximum allowable error 54>max — 10^^ mas, but also 
because they have periods very long, so that there is no hope they average out 
on reasonable Tofcs- Concerning the 18.6-year tide, by assuming an uncertainty of 
1.5% on fc2Q^^'^, the mismodeling on its perturbation amounts to -664 mas which is, 
however, very far from 5(j)„iax- 

4. Static geopotential perturbations 

As can be found in^^ , the perturbing function of degree I and order m of the static 
part of the geopotential can be cast into the form: 



^Ira — /II Flmp Glpq Sl/mpqj (^^) 



where: 

• Slmpq 



Clm 
— Sir, 



-ni even 



COS ^Impq 

-m odd 



-m even 



sin Iplmpq 

-m odd 



Sin 
Clr 

Cirri, Sim the unnormalized Stokes' geopotential coefhcients^^ 
• ipimpq — {I — 2p)uj + {I — 2p + q)M + m{Q — 9) in which 6 is the sidereal angle 

Concerning the preliminary analysis of the even degree perturbations, the same 
conclusions of the previous sections hold. All the perturbations on the radial direc- 
tions vanish, Arstatic = 0, while for the satellite' s azimuthal location only the zonal 
contributions are to be considered. Let us work out explicitly the perturbation due 
to the main even zonal coefhcient 6*2,0 • For the precessional secular rates induced 
by it of a satellite on circular orbit we have: 

dw 3 nC2.oR(B /I r ■2n /o1^ 
- = -^^(l-5cosz), (21) 
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? = ^^%^cosz, (22) 
at 2 

— = --— ^(3cosz -1). (23) 
By inserting eqs.(pl|)-(p3|), evaluated for i=0, into eg. (pi]) we obtain: 

^ = -■i./GM^C2.,oRla-^'^ = (2.637 • lO^^ cm^/^s-i) x a'^/^ (24) 

For a circular orbit of radius 7000 km the azimuthal secular rate is: 

^ = 1.89- 10^° mas/y. (25) 
at 

It is important to evaluate the error induced on such rate by the poor knowledge 
of the Earth' s gravitational field. According to the EGM96 model^^, a relative 
uncertainty of 7.3- 10~* weighs on the C2.0 coefficient. This yields 5<t) ~ 1380 mas/y, 
which is equivalent to 54) ~ 2.5-10^^ mas per revolution, being Porfe = ir ~ 5.82-10'^ 
s for a = 7000 km. Such error is 1 order of magnitude greater than S(f)max — 10~^ 
mas. 

Concerning the perturbations of odd degree, let us consider in detail the most 
important geopotential harmonic of degree / = 3. By reasoning as in the previous 
sections, it results that the satellite' s azimuthal location is not perturbed by the 
I = 3 part of the geopotential spectrum. Concerning Ar, the perturbation corre- 
sponding to the combination / = 3, m — 1, p ^ 1, q = —1 does not vanish. The 
full expression for the geopotential perturbation on the eccentricity is given by: 

Aeimp, = ^ Fimp [-{I - 2p)] Simp,. (26) 

na^ipimpq a e 

In our case it yields: 

Arstat^c = ^ ^GM^Rla-'/^ x ^ x [^3,1 cos ^-311-1 + ^3,1 sin tAsii-i], (27) 
with: 

Tp3ii-i = + (28) 

^pert 

With a = 7000 km we have for the period and the amplitudes of the harmonic 
terms: 

Ppert = 8.81 • 10^ s, (29) 



AAr = (-8.481 • 10^ cm) x 



(30) 



^ '5'3,l 

According to^^ the mismodeling weighing on the Stokes' coefficients of interest 
amounts to 6C3,i = 1.5 • 10~^°, 683^1 = 1.3 • 10^^". This leads to a mismodeled 
radial perturbation: 

Svstatic — (1.3 cm) cos t/ish-i + (1.15 cm) sin V'311-1. (31) 
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It results to be 2 orders of magnitude greater than the allowable Srmax — 10^^ mm. 
However, it must be pointed out that such a mismodeled perturbation averages out 
on an observational time span Tobs which is an integer multiple of 1 day, its period 
amounting to 24.48 hr. The same conclusions can be drawn for the other higher 
odd degree nonvanishing radial perturbations. 

5. Discussion and conclusions 

In this paper we have explicitly calculated, by averaging over one orbital revolution, 
the most relevant perturbations Ar and due to the dynamical and static part 
of the Earth's gravitational field on the radial and azimuthal locations of a satellite 
placed in an equatorial, circular orbit with radius of 7000 km. Furthermore, we have 
compared the mismodeling induced on such perturbations by the poor knowledge 
of the parameters of the Earth' s gravitational field to the maximum errors per 
revolution Sr„iax — 10~^ mm and S(j)„iax — 10^^ mas allowable in order to detect 
successfully the gravitomagnetic clock effect. 

Concerning the radial direction, it is affected by harmonic perturbations induced 
by the odd degree part, both static and dynamical, of the Earth's gravitational 
field. If, from one hand, the related mismodeling is 2 orders of magnitude greater 
than Srmax, from the other hand it must be pointed out that such mismodeled 
perturbations average out on not too long time spans since their periods range from 
1 to 50 days. 

The situation for the azimuthal angle is different. It is acted upon by the even 
degree zonal harmonics of the Earth' s gravitational field. The zonal tides are very 
insidious since they induce perturbations with great amplitudes acting on very long 
periods, so that it is necessary to wait for several years in order to average out their 
mismodeled effect which are up to 4 orders of magnitude greater than 64>max- The 
1 — 2, TO = part of the geopotential induces also a secular drift on the azimuthal 
satellite' s angle; the uncertainties in 6*2,0 induces on it a mismodeled rate per 
revolution which is 1 order of magnitude grater than 64>max- 

The conclusions outlined here hold for r = 7000 km; let us see how the situation 
changes with different values for the orbital radius. The possible scenarios turn out 
to be very intricate. Indeed, from one hand we have secular or semisecular mismod- 
eled perturbations which could be reduced only by enlarging the orbit radius, and 
from the other hand there are the periodic mismodeled perturbations whose periods 
grow with the orbit radius making so much more difficult to average them out on 
reasonable time spans. Moreover, it should also be considered that the maximum al- 
lowable errors depend on the orbit radius and they decrease with increasing radius 
putting, in this way, more stringent constraints on the mismodeled gravitational 
perturbations. This is shown in Fig.(|l|). In Fig.(||) and Fig.(|^) we show how the 
mismodeled perturbative amplitudes due to the 18.6-year tide and the C2,o depend 
on the orbit radius. From an inspection of Fig.(|l|) and Fig.(||) it can be noted that 
the major problems come from the azimuthal error and the perturbation induced 
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Gravitomagnetic clock effect error budget 



'0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Satellite orbit radius (km) ^ 



_ 0.025 



T 0.02 - 
o 

n 0.015 - 
E 

I 0.005 - 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Satellite orbit radius (km) ^ ^ 

Fig. 1. Maximum allowable errors in the radial and azimuthal locations. The values for the orbit 
radius span from 7000 km to 42160 km for a geostationary satellite. 
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Fig. 2. Mismodeied azimuthal perturbation induced by the 18.6-year tide. 
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Secular mismodeled azimuthal perturbation per revolution 
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Fig. 3. Mismodeled azimuthal rate per revolution induced by the C2,o geopotential coefficient. 
The values for the orbit radius span from 7000 km to 42160 km for a geostationary satellite. 

by the 18.6-year tide: as the orbit radius grows, the mismodeled tidal perturbation 
is always greater than the maximum allowable error by 2 or 3 orders of magnitude. 
Since this important source of systematic error cannot be made harmless by vary- 
ing the orbital parameters of the satellites, it should be necessary to average out its 
effect: but this means to choose a time span Tots — 18.6 years at least. 

By inspecting Fig.(||) it can be noted the growth of the period of the radial tidal 
perturbation is induced by the Ki, I = 3 p ^ I q = —I. This is an important 
feature since its mismodeled amplitude is at cm level and is independent from the 
orbit radius. Moreover, Sr^ax is of the order of 10~^ mm, so that we could eliminate 
the effect of such a perturbation only by averaging it over an integer multiple of its 
period. 

We can conclude that the present level of accuracy in the knowledge both of 
the Earth solid and ocean tides, and of the static part of the geopotential docs not 
allow an easy detection of the gravitomagnetic clock effect at least by using short 
arcs only. 
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X 10"* tidal perturbation period 




Satellite orbit radius (km) ^ ^ q* 

Fig. 4. Period of the Ki I = 3 p = 1 q = —1 radial tidal perurbation. The values for the orbit 
radius span from 7000 km to 42160 km for a geostationary satellite. 
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